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ON ELECTROMAGNETIC SWITCHING PROPERTIES OF A BI-STABLE FLUIDIC ELEMENT 

WITH NON-DESTRUCTIVE MEMORY 

By Joseph S. Koziol, Jr., and Robert De Furia 
Electronics Research Center 


SUMMARY 

This report discusses the feasibility of controlling a 
magnetically permeable liquid bead with an electromagnet between 
two stable states of a particular system. Several useful proper- 
ties of the switching system, namely, magnetomotive force re- 
quired for switching, switching time, and ambient acceleration 
field limits are determined in terms of the system parameters. 

This concept suggests a simple, bi-stable, low-power, 
el ectri c-to-f 1 uid transducer with non-descructi ve memory. (The 
fluid bead retains its last stable state with all power off.) 

Development of a magnetically permeable liquid with high 
surface tension did not reach the stage for implementation, at 
the completion of the analysis, to permit experimental verifica- 
tion of the analytical results contained herein. 


INTRODUCTION 

The unique characteristics of fluidics offer many advantages 
over conventional methods for performing logic, sensing, amplify- 
ing, and actuating. Besides eliminating wear of moving parts, 
fluidic implementation is relatively impervious to damage by the 
environmental effects of radiation and high temperature. These 
considerations are especially important when designing systems 
for deep-space missions such as a space probe. 

Recognition of these inherent advantages, per se, has not 
and will not be a panacea for realizing actual functioning sys- 
tems. Although fluidic technology is expanding at a rapid pace, 
system development is still in its embryonic stage and many break- 
throughs are needed. One such area of significance is the effi- 
cient conversion of electric signals to fluidic equivalents. 

The purpose of this paper is to describe a system with this 
capability and to discuss the feasibility of its implementation. 

The mechanism analyzed herein consists of a liquid bead 
which is geometrically retained in either of two stable states by 
its surface tension. In a typical system, each state could serve 
as a control port for gas flow. 


This concept of controlling a gas flow by means of a liquid 
bead has been labelled, perhaps inappropriately, two-phase fluid- 
ics. It was developed by the Conrac Corporation and now provides 
a basis for research on a number of fluidic components. Its use- 
fulness has already been demonstrated by the successful operation * 
of an alpha-numeric display with memory (refs. 1 and 3) and by an 
infinite-input impedance amplifier (ref. 2). In these two appli- 
cations, control of the liquid bead was accomplished by pressure , 
signals. This report considers control of a magnetically perme- 
able liquid bead by means of an electromagnet; hence, suggesting 
possible applications as an electric to fluid transducer. 

This report evolved from previous research originated by 
the Conrac Corporation (formerly Giannini Controls Corporation) 
under Contract NAS 12-43 with NASA-ERC (refs. 1 and 2). 


SYMBOLS 

O 

A cross-section area of electromagnet, meter 

2 

Ag surface area of bead, meter 

2 

A q initial surface area of bead, meter 

2 

B magnetic field intensity, weber/meter 

D height of electromagnet gap, meter 

Fg electromagnetic force, (amp-turn) weber/meter 

Fgy surface force, newton 

2 

g earth's gravitational constant, meter/sec 

g. ambient acceleration field factor, dimensionless 

a 

H bead dimension, meter 

I current, ampere 

kg kilogram 

£ electromagnet length, meter 

L restriction length, meter 

M mass of bead, kilogram 

N number of turns of coil 


2 


NI ,MMF 


magnetomotive force, amp-turn 
radius of restriction, meter 
radius of chamber, meter 


r 

R 


R > R g' R l 
R 2 > R 4 

t 


T 

U 

U, 


’BO 

*ST 

J^STO ’ 

J ST0 


W 

X 


x sw 

* 


x 


P 

* 

a 


reluctance, ampere/weber 

time, second 
switching time, second 
kinetic energy, newton-meter 
potential energy, newton-meter 

electromagnetic potential energy, (amp-turn) weber 
electromagnetic energy constant, (amp-turn) weber 
surface potential energy, newton-meter 

surface energy constant, newton-meter 

electromagnetic energy in gap, (amp-turn) weber 
electromagnet width, meter 

independent coordinate, center-of-mass position of 
bead, meter 

position of bead's leading surface, meter 
geometric center of switching system, meter 
position of center-of-mass at switching, meter 
velocity of bead's center-of-mass, meter/sec 

2 

acceleration of bead's center-of-mass, meter/sec 

3 

density of bead, kg/meter 
magnetic flux, weber 

surface tension constant, newton/meter 



magnetic permeability of electromagnet, weber/amp-meter 

magnetic permeability of bead, weber/amp-meter 

magnetic permeability of body in electromagnet gap, 
weber/amp-meter 

magnetic permeability of air, weber/amp-meter 


ANALYSIS 


The generalized Lagrange equation is a convenient method 
for analyzing the dynamic behavior of discrete point systems 
* especially when the system can be described in terms of independ- 
ent coordinates. It will be used to describe the motion of a 
magnetically permeable liquid bead constrained to move through a 
, restriction while influenced by a magnetic field. The general 
representation of such a system is shown in Figure 1. 

In analyzing the bead's motion, emphasis will be placed on 
the center-of-mass of the bead, while the deformation of the bead 
will be assumed to affect only the potential energy. In other 
words, the bead is treated as a rigid body with its kinetic energy 
balanced by the electromagnetic energy of the applied "B" field 
and surface energy, or work done by the surface tension, which is 
implemented by the variable restriction. 

A strict energy description of the problem would require 
representing the variational motions of each individual segment 
of the fluid bead with properly defined coordinates and constraints. 



Figure 1. -General representation of a magnetically 
permeable bead in a variable restriction 
with applied magnetic field "B" 


However, since the main point of interest is not to examine the 
fluid flow phenomena involved, but rather to follow the macro- 
scopic displacement of the bead, the net effect will be repre- 
sented in terms of its center-of-mass motion. 
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This approach allows the problem to be described in terms 
of a single independent coordinate, namely, the motion of the 
bead center-of-mass . 

The Lagrange equation for one independent coordinate, x, 
is given by: 


_d_ ( . 3(T - U) = o 
dt \ 3 X/ 3 X 


( 1 ) 


where T and l) are the total kinetic and potential energies of the 
system, respectively. 

No generalized force is included on the right-hand side 
of this equation, since the applied electromagnetic force arises 
from a potential. The potential energy U , thus, consists of the 
surface tension energy and the electromagnetic energy. 

In order to proceed further in the analysis, a specific 
system needs to be defined. The system chosen for this report, 
shown in Figure 2, lends itself to a simple description of the 
restriction parameters and is not intended to suggest an optimum 
design . 

The origin of the coordinate system is established at the 
bead's initial center-of-mass position when the bead is just 
about to enter the restriction, while "x 0 " is the distance to the 
system's geometric center from the origin. 



Figure 2. -System configuration 
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An electromagnetic energy approximation and a surface 
tension energy approximation have been developed in Appendices A 
and B, respectively, for a certain switching model of the system 
shown in Figure 2. These approximations are repeated here for 
* convenience. The electromagnetic energy approximation is: 




where 


U B U B0 


+ FgX 


( 2 ) 


’BO 


4(ni) 2 < 


[ y o( 1 'ir)( 2R+L )] + [ M o iH 2R+l )] 



(3) 



P b and y Q are the magnetic permeabilities of the bead and air, 
respecti vely . 
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The surface energy approximation is: 


U ST U ST0 + 

f st x 

X < X 

^STO 

f st x 

X > X 


where 

U^jq = 2ttRo(R + H) 

IL Tn = energy constant determinable 
from Figure B-4 



( 5 ) * 


( 6 ) 



and 

H + I 

x Q = — 2 — = geometric center. 

The kinetic energy of the bead is: 

T = 1 Mx 2 


(7) 


( 8 ) 


(9) 


where M is the mass of the bead and x is the velocity of the 
bead's center-of-mass . 

When these equations are substituted into Eq. (1), the 
equation describing the behavior of the bead becomes 


Mx = -Fg - F $t 

X < X 

" f b + f st 

X > X 


( 10 ) 
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The value of -Fg is positive, since Fg is negative and 
represents the electromagnetic driving force contribution result- 
ing from an applied magnetomotive force, NI. 

The quantity -F ST represents a constant restraining force 
as the bead enters the restriction (x < x 0 ) and Fcf represents 
an equal propelling force when it emerges (x > x 0 ). 

The solution of Eq. (10) describes a single switching 
process and assumes that friction is negligible. Further, the 
switching process is treated as symmetrical and reversible; thus, 
impact effects are neglected. 

The solution of Eq. (10) is valid only when the right-hand 
side of the first equation is positive. Otherwise, the applied 
electromagnetic force, -Fg, is less than the surface tension 
restraining force and the bead will not switch. In this case, 
the bead has no motion and any applied electromagnetic force is 
balanced by a wal 1 -restrai ni ng force. 

When the applied electromagnetic force is equal to the sur- 
face tension restraining force, the bead is on the verge of 
switching, and this defines the electromagnetic force required 
for switching. 

If the applied electromagnetic force is greater than the 
electromagnetic force required for switching, then switching time 
can be determined from solution of Eq. (10). 

The determination of the ambient acceleration field that 
will cause switching without an applied electromagnetic force 
will provide a third characteristic of the system. 

These three characteristics, when related to the system 
parameters, would provide valuable design information. They will 
now be examined separately. 


Electromagnetic Force Required for Switching 

The electromagnetic force required for switching can be 
determined, by definition, from the following equation. 

-Fg = F $t (11) 


Substituting from Eqs. (4) and (7) and rearranging the 
equation in normalized form gives: 
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02 ) 



where NI is the magnetomotive force. 

This curve is plotted in Figure 3 as a function of the 
restriction r/R. It can be seen that the magnetomotive force 
required for switching is independent of the size of the bead, 
H/R, the restriction length, L/R, and the permeability ratio, 
u b /u 0 y b / y o > 1 0 0 ) • 

As would be expected, the electromagnetic force required 
for switching decreases with larger restrictions. The curve also 
shows that for blocked restrictions ( i . e . , r/R = 0), it is impos- 
sible to switch the bead, and for no restriction (i.e,, r/R = 1), 
the bead can switch without an applied magnetomotive force. 


Switching Time 

It is first necessary to define "switching." "Switching" 
will be defined to have occurred when the trailing surface (see 
Figure 2) reaches the restriction. Switching time can now be 
determined from the dynamic solution of Eq. (10). 

Substituting from Eqs. (4) and (7), solving, and normali- 
zing properly gives: 



+ 


R 

r 




(contd) 
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Plots of Eq. (13) are shown in Figures 4, 5, and 6 as 
functions of the applied magnetomotive force. The vertical 
asymptotes are the magnetomoti ve forces required for switching 
for each respective restriction. 


1 1 



JJ b /iJ 0 -100, 200. 1000 



iotive force required for 
ig vs. restriction size 


r/R=0.5 H/RHO 


L/R *1.0 



Figure 5. -Switching time vs. applied magnetomotive force 
showing the effect of the bead's size 



Figure 6. -Switching time vs. applied magnetomotive force 
showing the effect of restriction size 
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As would be expected, for a given restriction, higher 
applied magnetomotive forces result in lower switching times. 

Figure 4 shows that the effect of restriction lengths is 
negligible for L/R<1. 

Figure 5 shows that the size of the bead can reduce the 
switching time significantly. 

Figure 6 shows that for a given applied magnetomotive force 
larger than that required for switching, the switching time de- 
creases significantly with larger restrictions. 


Ambient Acceleration Field Limit 

It is possible that ambient acceleration fields could cause 
unwanted switchings without an applied magnetomotive force. The 
axial acceleration field component that would cause such switch- 
ing can be determined from Eq. (10) with the magnetomotive force 
set equal to zero. 

Substituting from Eq. (7) and normalizing properly gives: 



Equation 14 is plotted in Figure 7 as a function of the 
restriction, r/R. 
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H/R* 



(bd^/^)/ D b 
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Figure 7. -Ambient acceleration field limit vs. restriction size 



It is apparent that the curves are not valid for r/R<.25. 
This can be accounted for by an increased error resulting from 
the surface tension energy approximation in this region. (See 
Append i x B . ) 

Otherwise, it can be seen that the axial acceleration field 
limit is independent of the restriction length and inversely 
related to the size of the bead. Also, as would be expected, the 
bead could switch with virtually no ambient acceleration field 
present when r/R = 1 . 


Numerical Example 

A numerical example is provided here to show that the pre- 
dicted values of magnetomotive force required for switching, 
switching time, and axial acceleration field limit are practical. 
Consider a typical system with the following parameters: 

L/R = 1 .0, r/R = 0.5 
S/R = 1.0 

y b /u 0 = 200 

= 4 7T x 10~ 7 weber/amp-meter 

_ "5 

R = 3.0 x 10 meter 
a = .462 newton/meter 
p = 13.5 x 10 3 kg/m 3 

Magnetomotive force required for switching .- From Figure 3: 



(NI) = 88 amp-turns. 


For N = 1000, and a resistance in the coil equal to 1 ohm, 
the power expended during switching is .0077 watt. 
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Switching time .- From Figure 5 and with an applied mag- 
netomotive force of 200 amp-turns: 


t 


s 


2R 



.94 


t = .051 second. 


If the applied magnetomotive force were increased to 500 
amp-turns, the switching time would be reduced to only .028 sec- 
ond . 


Ambient acceleration field limit .- From Figure 7: 



g a = .44 g ' s. 


CONCLUSIONS 

Several useful properties of the electromagnetic switching 
system have been determined in terms of the system parameters. 
These properties are magnetomotive force required for switching, 
switching time, and ambient acceleration field limit. 

The results of the numerical example predict that the 
system could perform quite satisfactorily as a low-power, digital, 
el ectri c-to-f 1 ui d transducer. 

It is quite difficult to otherwise appraise this system 
without experimental data. Several assumptions were made in the 
development and could be validated by actual experimentation. 
However, the development of a proper fluid bead has not yet 
reached the stage for final implementation of this concept. 
Experiments with various materials for developing a non-wetting 
liquid with high magnetic permeability and high surface tension 
are progressing in this Center. 
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The results of this report are not valid for small restric- 
tions, first of all, because of the approximation method adopted. 
(See Appendix B.) Also, the actual switching model equation omits 
the effects of friction and hysteresis which would be significant 
for small restrictions. 

No attempt was made to perform an accurate small restric- 
tion analysis because such a system appears to be impractical. 

In addition, the analysis would require a much more complicated 
model to describe the behavior of the bead which would be subject 
to greater possibilities of rupture. 

Regardless of these comments, the fluid bead concept itself 
might suggest other applications and thereby warrants further 
research to expand fluidic technology. 

The fabrication of a proper liquid bead is a primary need. 

The required properties include high magnetic permeability, non- 
wetting capabilities, and high surface tension. True fluids 
having high magnetic permeability have been developed but these 
are wetting fluids and, therefore, unsuitable for this applica- 
tion (refs. 4 and 5). In most cases, they are corrosive, highly 
volatile, and prone to oxidation. Also, surface tension is often 
treated as a property of liquids, though it is actually a property 
of the boundary of the liquid and cannot be given an exact value 
without specifying the material which lies on the other side of 
the boundary. Hence, development of the bead material cannot be 
separated from that of the system material. 

Other germane research areas could include optimization of 
the restriction geometry to allow switching with minimum power 
consumption and optimization of performance (i.e., quickest 
switching without damage) by application of bang-bang control to 
the applied magnetomotive force. 
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APPENDIX A 


DERIVATION OF ELECTROMAGNETIC ENERGY 
MODEL AND APPROXIMATION 


In simple magnetic circuits, as shown in Figure A.l, there 
is a formal analogy to Ohm's law for electric circuits given by: 

MMF = R<t> (A.l) 


where MMF is the magnetomotive force analogous to electromotive 
force, R is the reluctance of the material in the circuit analogous 
to resistance, and <j> is the flux analogous to current. 

The magnetomotive force may be regarded as the agent that 
establishes the flux. The source of the magnetomotive force is 
the current flowing in the coil wound around the core of the mag- 
net and is given by: 

MMF = N I (A. 2) 


where N is the number of turns of the coil and I is the current in 
the coil. The dynamics of the electronic circuit are neglected. 



Figure A . 1 . -El ectromagnet 

Reluctance is defined as the ratio of magnetomotive force 
to the flux it produces. For a straight bar of uniform cross- 
section A, and length D, through which the flux lines are straight, 
parallel, and uniformly distributed (no fringing), the reluctance 
is given by: 
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I 


R = 


D 

yA 


(A. 3) 


in which u is the magnetic permeability of the material. 

Since the permeability of the ferromagnetic core is much 
larger than the air or vacuum gap, the total reluctance of the 
magnetic circuit can be reasonably approximated by: 



(A. 4) 


in which yg is the permeability of the material in the gap, 
as long as D is small . 

Combining the above equations: 


NI = R g <J> . (A. 5) 


The energy in the gap can also be expressed in terms of the 
reluctance of the magnetic path. With non-linearity and hysteresis 
neglected, the energy is: 



(A. 6) 


Now if a body of permeability yg is placed in the gap, as 
shown in Figure A. 2, the equivalent reluctance of the circuit 
and the energy in the gap will be affected. 


The equivalent reluctance in the gap analogous to electronic 
circuitry becomes that of two rel uctances , D/ (y Wx) and D/[y W(a-x)] 
acting in parallel, namely: 9 0 


v = JL. 

K g £W 


.-gif)' w o 0 -f) 


(A. 7) 


where the area in the gap is now assumed rectangular and W is the 
width and i is the length. 
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K 


Figure A . 2 . -El ectromagnet showing body (with 
permeability y^) partially in gap 

This is also equivalent to the reluctance of the entire 
circuit if: 


u » u. 


(A. 8) 


From Eq. (A. 6) the energy in the gap is: 


u g ' 4 * 2 *g - 4< NI >7 


(A. 9) 


Og = -^(NI ) 2 ^[w g X + « 0 <^-x)] • 


(A. 10) 


A similar process can be followed in determining the energy 
in the restriction of the system shown in Figure A. 3. 


Lk- 




El ectro- 
magnet 



Figure A . 3 . -Comp! ete switching system 




In this configuration, the equivalent reluctance consists 
of parallel and series combinations schematically represented in 
Fi gure A . 4 . 



Figure A . 4 . -Schemati c of Electromagnetic Circuit 


The reluctances are defined as follows: 



where and y 0 are the magnetic permeabilities of the bead and 
air, respectively, xl is the position of the leading surface, and 
L is the restriction length. 

These equations assume that the permeability of the material 
containing the bead is approximately equal to that in the empty 
cavity (air) and that the bead has a flat leading surface. 

Only one electromagnet is shown in Figure A. 3, whereas 
in the actual system a second electromagnet, for reverse switch- 
ings, would be arranged similarly on the left side but rotated 
90 degrees about the switching axis to prevent mechanical inter- 
ference. 
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When the equivalent reluctance is determined and combined 
with Eq. (A. 9), the energy in the restriction becomes: 


+ [ u b tJ o rx L ~ |1 

ho + u b (R * r) Jj 


(A. 15) 


Equation (A. 15), therefore, relates the electromagnetic 
potential energy of the bead to the position of its leading sur- 
face X|_, as the bead proceeds through the restriction. The form 
of Lagrange's equation in the text requires that the potential 
energy be expressed as a function of the independent coordinate, 
namely, the position of the bead's center-of-mass . 

The position of the bead's leading surface can be related 
to its center-of-mass position by imposing the constant volume 
constraint. This gives: 



where xl and x are defined from the origin of the coordinate 
system as shown in Figure A. 5. The origin of the coordinate 
system is established at the bead's initial center-of-mass posi- 
tion when the bead is just about to enter the restriction. 

Now consider the bead as it exits the restriction. If the 
switching process resembled the "squirt" model shown in Figure 
A. 6, the relationship between the position of the bead's leading 
surface and its center-of-mass position would still be described 
by Eq. (A. 16). However, deviations from the "squirt" switching 
model lead to different x^ vs. x relationships. 
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Trailing Surface 


Center-of-Mass 


Figure A. 5. -Bead configuration showing leading surface, 
trailing surface, and center-of -mass 




The bounds of any possible vs x relationship can be 
represented quite simply with respect to the "squirt" model. 
Equation (A. 16) is plotted in Figure A. 7 as solid lines for 
several restriction sizes. The leading surface position is the 
ordinate and the center-of-mass position is the abscissa. 

Note that each curve can be approximated by a straight 
line quite well as x/R becomes large. 

Now for each respective restriction, the solid line repre- 
sents the maximum values of x^ that are physically possible. 

Conversely, the straight line of slope "1" would represent 
the minimum values of xl that are physically possible for any 
restriction (i.e., the center-of-mass and leading surface have no 
relative motion). The straight line of slope "1" also represents 
the only possible x^ vs. x relationship when r/R = 1. 
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Figure A . 7 . -Leadi ng surface position vs. center-of-mass 
position for several restriction sizes 



The former assumes that the bead does not contract more 
than the restriction size while switching and the latter assumes 
i ncompressi bi 1 i ty . 

Accordingly, it can be deduced from Figure A. 7 that: 

(1) Equation (A. 16) is a good, reasonable description of any 
switching process for large restriction sizes (i.e., r/R > .75). 

(2) The dotted straight lines shown are good approxima- 
tions to the solid lines for x < x Q (geometric center). 

Therefore, rather than involve the non-linearities of and 
undetermined deviations from the "squirt" switching model, the 
straight-line approximations are used in the text in analyzing 
the system characteristics. 

These approximations should be extremely accurate for pre- 
dicting the ambient acceleration field limit and the magnetomotive 
force required for switching, since only the initial x L vs. x 
information is required. That is, these properties depend only 
on the largest rate of energy change and this occurs in the 
restri cti on . 

The approximations are not so accurate for determining the 
switching time since significant variations occur before complete 
switching. Also, the approximations are not conservative (they 
tend to predict switching times less than the actual switching 
times), but still predict a reasonable estimate (i.e., within a 
factor of 2 for r/R = .5). 

The straight-line approximation relating the bead's leading 
surface position to its center-of-mass position can be determined 
from Eq. (A. 16) and is given by 


x 


L 


H 

2 


+ 





2 



(A. 17) 
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Finally, the electromagnetic energy can be determined as a 
function of the bead's center-of-mass position by substituting 
Eq . ( A . 1 7) into Eq . ( A . 1 5 ) : 

U B = U B0 + F b x (A. 18) 


where 


and 


'BO 


4< NI > 2 {to( 1 -B-)( 2R+L )] + [ u oiH 2R+L )] 




(A. 19) 



+ 






(A. 20) 


The electromagnetic energy is also affected by possible 
deviations in reluctance as the bead emerges from the restriction. 
But, this effect tends to offset the x L vs. x deviations just 
mentioned. This can be shown by examining Eq. (A. 15). Any devia- 
tion in reluctance would mean an effective increase in r and, 
hence, an increase in the stored electromagnetic energy, which, 
in turn, would reduce the switching time. 
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This effect can be examined quantitatively for the hypo- 
thetical switching process model shown in Figure A. 8. This model 
reflects the minimum possible reluctance across the chamber as a 
function of the bead's leading surface. 


// /////. 



Center-of-Mass 


s//'7 /~// // 


Figure A . 8 . -Hypotheti cal switching process model 


After computing the effective reluctance in the restriction 
and right chamber as a function of the bead's leading surface, 
and placing into the energy Eq. (A. 9), the electromagnetic energy 
becomes : 


U 


B 



+ 


u o y b 


M o r + 


rL 


% 


L 


+ p Q ( 4 R + L ) + 



(A. 21) 


In order to asses the effect of reluctance on switching 
time, the coefficients of X|_ in Eqs. (A. 21) and (A. 15), which 
represent the energy rate of change, must be compared. 

It is found that the coefficient of the hypothetical switch- 
ing process is larger than that of the "squirt" switching process 
model by a factor of 
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According to the switching time equation in the text, this 
reduces the estimated switching time by a factor of approxi- 
mately: 


1 r y r\r ) 


This factor is shown in Table A.l for different values of 
restriction sizes and permeability ratios. 

This development was intended to show that the estimated 
switching time given by the "squirt" switching process model 
tends to be conservative and not intended to determine a minimum 
switching time. The switching time predicted by the "squirt" 
switching process model should provide a reasonable estimate for 
any practical switching system. 


TABLE A.l 

ELECTROMAGNETIC SWITCHING TIME FACTOR 


r/R 

u b /u 0 = 1f)0 

u b /u o = 200 

1 .0 

1 

1 • 

.50 

14 

20 

.33 

25 

35 

.25 

35 

49 
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APPENOIX B 


DERIVATION OF SURFACE TENSION ENERGY 
MODEL AND APPROXIMATION 


The energy expression for surface tension is given by 

U ST / ad Ag + constant (B . 1 ) 


where U$y is the potential surface energy, Ag is the surface 
area of the bead, and a is the surface tension constant. 

It may appear, at first, that the expression relating the 
surface area of the bead to its center-of -mass position is simply 
a geometric problem. A more thorough examination will reveal 
that the bead can exit the restriction with various orofiles. 

In this analysis, an energy model will be depicted for a 
certain type of switching process. It will then be shown that 
any practical switching process can be described reasonably well 
by this model. Finally, a linear approximation will be fitted 
to the model to allow simple yet credulous analysis of the 
system characteristics. 

Consider the bead shown in Figure B.l. The menisci of the 
bead are assumed flat to permit an uncomplicated description of 
the bead's surface area. Furthermore, the thermodynamic crite- 
rion of marginal stability states that variations in the meniscus 
shape produce no variations in the total potential energy (ref. 
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The origin of the coordinate system is established at the 
bead's initial center-of-mass position when the bead is just 
about to enter the restriction. 

The relationship between the surface area and the center- 
of-mass position of the bead as it proceeds through the restric- 
tion, as shown in Figure B.2, is given by: 


A b = 2irR ( R+H ) 




+ 



+ 



(B.2) 


This relationship is derived by first relating the surface 
area to the leading and trailing surfaces of the bead; relating 
the center-of-mass position to the leading and trailing surfaces; 
imposing the constant volume constraint; and finally, combining 
these expressions. 


New Center-of-Mass 
Position 



Figure B.2. -Bead proceeding through restriction 
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Now consider the bead as it exits the restriction. If the 
bead were forced strongly to switch, the switching process would 
resemble the "squirt" model shown in Figure B.3. Only during 
this type of switching process could the bead attain its maximum 
« energy (area) state. Obviously the area relationship of Eq. (B.2) 
still applies here. 



Now, rather than trying to depict deviations from the 
"squirt" switching process physically, it is more convenient and 
apropos to examine deviations from a plot of the "squirt" switch- 
ing process [Eq. (B.2)], which is shown as the solid line in 
Figure B.4. The surface area of the bead, A B , normalized with 
respect to its initial area, A 0 » is the ordinate, and the center- 
of-mass position, normalized with respect to R, is the abscissa. 

Because of the symmetry of the system, the complete single 
switching process is shown symmetrical about the geometric center. 

The solid line describes any switching process as the 
leading surface of the bead proceeds through the restriction. 

Only when the bead begins to emerge from the restriction can 
deviations occur. The shaded portion for each respective re- 
striction represents the possible deviations from Eq. (B.l) for 
that restriction, as the bead emerges from the restriction. 

In effect, the minimum value of each shaded portion repre- 
sents the minimum possible surface area (energy) of the bead at 
switching for the particular restriction. It is the area of the 
bead when the leading surface just begins to emerge from the 
restriction, as shown in Figure B.5. 
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Figure B . 5 . -Representati on of the minimum possible area 
that the bead can have at switching 


It is apparent from Figure B.4 that: 

(1) Equation B.2 is a good reasonable description of any 
switching process for large restrictions (i.e., r/R>.75). 

(2) The dotted straight lines shown are good approxima- 
tions to the solid lines. 

Hence, rather than involve the non-linearities of and 
undetermined deviations from Eq. (B.2), the straight-line approxi 
mations are used in the text to analyze the system characteris- 
tics. These approximations are especially warranted in deter- 
mining the ambient acceleration field limit, and magnetomotive 
force required for switching, since these properties depend only 
on the largest rate of energy change and this occurs in the 
restriction. 

While not entirely as accurate for determining the switch- 
ing time, the straight-line approximation does tend to give a 
reasonable, conservative estimate. 

Using Eqs. (B.2) and (B.l), the straight-line approximation 
relating the surface energy of the bead to its center-of -mass 
position becomes: 


= U 


where 


U< 

U, 


STO + F ST x 

x < x„ 
0 

( B . 3 ) 

STO " F ST X 

x > x„ 
0 


2uRa(R+H) 


(B.4) 

energy constant 

determi nabl e 

from Figure B.4 
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r 



+ 



1 



( B . 5 ) 


and 


x 


o 


H + L 
2 


geometric center. 


( B . 6 ) 


It is interesting to note that the "squirt" switching model 
curves begin to decrease for large x/R, as r/R becomes less than 
.25. This effect is shown in Figure B.6. What this means is 
that the center-of-mass of the bead is now more affected by the 
large moment arm of the mass past the restriction than by the 
mass remaining behind in the chamber. 

Also, note that the initial slopes of all the switching 
model curves continue to increase for decreasing values of r/R. 
These slopes determine the magnetomotive force required for 
switching and the ambient acceleration field limit and, to a great 
extent, the switching time. The monatonic nature of these slopes 
suggests that the "squirt" model is valid for all restr i cti ons . 
However, since the initial slope information is lost completely 
in the approximation, it must be concluded that the properties of 
the switching system determined in the text are not valid for 
r/R <.25. 

It remains, though, that the switching model curves are 
theoretically correct for all restrictions under the assumptions 
stated (i.e., no friction, hysteresis, etc.). 
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Figure B. 6. -Surface area vs. center-of-mass position showing effect 
of small restriction sizes 
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